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It is well known that the grand unification group £77(5) relates the mass matrices of the 
down quarks and charged leptons. There is some empirical support for the existence of such 
a relationship in the fact that when the fermion masses are extrapolated to the GUT scale 
in the MSSM one finds mj = m T WW j m s = m M /3, and rrid = 3m e j3|. However, the 
pattern of masses of the up quarks is very different. One difference is that the t mass is 
much greater than the b and r masses, which is usually explained by saying that the ratio 
of VEVs v u /vd = tan/? is large compared to one. Another difference is that the interfamily 
mass hierarchies are much stronger for the up quarks than for the down quarks and charged 
leptons (e.g. m c /m t <C m s /mb and m u /m c <C m^/ms). It is tempting to say that the up 
quark mass matrix (Mu) is more distantly related to the down quark and charged lepton 
mass matrices (Mjj and Ml) than the latter are to each other. On the other hand, there is 
the tantalizing fact that m u /md ~ 1. 

In this paper we suggest a somewhat new approach which qualitatively explains why 
m t ^> m& but m u ~ m^. The idea is that there are "underlying" mass matrices (denoted 
by the superscript zero) whose structure is controlled by £0(10) and which satisfy M v ~ 
~ M° ~ (it is assumed v u /vd ~ 1), but that a strong mixing of the third family 
with vectorlike fermions at the GUT scale distorts these underlying mass matrices in such 
a way that m& and m T are highly suppressed relative to m t . This distortion does not affect 
the first family much, so the masses m u , rrid, and m e remain of the same order. 

The approach we will describe has several other virtues: (a) It can be realized in models 
with very few parameters, (b) It dovetails with the ideas of Ref. |4J for solving the Strong 
CP Problem. And (c) it implements the "lopsided" mass matrix approach to explaining 
large neutrino mixing angles j^, 0| . 

To understand the idea, consider an £77(5) model (which will later be assumed to descend 
from £0(10)) that has three families of fermions in 10j + 5j + lj, with mass terms of the 
form (M^ijlOilOj + (M^) ii 10 i 5 i + (M£) ii 5 i 10 j + (M^yM,-. The matrices M% and M° L 
come from the VEV of a 5 of Higgs, and My and M% (the neutrino Dirac mass matrix) come 
from the VEV of a 5 of Higgs. Suppose that there are also for each family a vectorlike pair 
of quark/lepton multiplets, denoted 5^ and 5' { and having superheavy mass terms AySfoj + 
Bij5i5'j. (Aij ~ Bij ~ Mqut-) There is then mixing between the ordinary three families 
and the vectorlike fermions, more specifically the mixing is between the 5; and the b[. (A 
similar idea, but with mixing among fermions in 10's of £77(5) was used in [5J]. However, 
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the models proposed there were very different in character from the present models.) With 
the mass terms specified above, we may write 



5 5^M F 



10 \ ( (Mp)ij B i3 









\ 5 i ) 



where F = L or D, and M£i is either M° or (M^>) T , depending on whether the fermions in 
5 are or g?£. 

In order to find the light fermion mass matrices in the effective low-energy theory, we 
must do a unitary transformation Aip = UAi° F that eliminates the off-diagonal block B in 
the full mass matrix given in Eq. (0). Such a transformation is 



W = „ T h ( 2 ) 




where x = A = (/ + xxt) -1 / 2 , and A = (/ + x^x)~ 1 ^ 2 . (To check the unitarity of Li 

it is useful to note that x^A = Ax* and a; A = Ax.) This gives the result for the low energy 
mass matrices 

M L = AMI (3a) 
M D = M° D A T . (3b) 

Basically, the hermitian matrix A describes the mixing of 5j with 5^ It appears on the left 
in the equation for Ml since (Ml)^ couples to 5jl0j . It appears on the right in the equation 
for Mr, since (Mc)ij couples to 10j5j. For the Dirac neutrino masses we have 

(«)-(;,)=(«)(«'• ;)(;;)• 

giving 

M N = AM". (5) 

Since the masses of the up- type quarks come from a lOjlOj coupling of the fermions, they 
are not affected by the mixing of the 5; with the 5j. Consequently, 

M v = M°. (6) 

Before we discuss how the structure we have described can help us explain the magnitudes 
of quark and lepton masses and mixings, we note that it is exactly the kind of structure that 



is used in the solution of the Strong CP Problem proposed in ^|. The idea there was the 
following. Suppose that CP is a symmetry of the lagrangian that is spontaneously broken, 
and that the VEV that breaks CP appears in the off-diagonal matrix B in Eq. JI}, but 
not elsewhere in the quark mass matrices. Then and A are real, and it is easily shown 
that the determinant of the full mass matrix A4d is therefore real. Also real, of course, 
is the determinant of My. Thus, at tree level, the phase 9 is zero. At higher order, these 
matrices can receive complex corrections that induce a non- vanishing 9, but these may be 
made small. (In SUSY, there can be contributions to the 9 parameter that are harder to 
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make small, for example, one-loop corrections to the gluino mass p|. How large these are 
depends upon how SUSY is broken. These contributions are not a problem in theories with 
gauge-mediated SUSY breaking, for example. We imagine that whatever mechanism resolves 
the usual SUSY flavor and SUSY CP problems will also suppress these extra contributions 
to 9.) On the other hand, since B is a complex matrix, so is the matrix x = BA^ 1 and 
the matrix A = (/ + xx^)~ 1 ^ 2 . Consequently, the mass matrix of the light three families of 
down-type quarks in the effective low-energy theory, given by Md = M^A T , is also complex, 
which means that in general there is a non- vanishing Kobayashi-Maskawa phase. 

In short, the structure in Eq. allows a spontaneously generated phase in the matrix 
B to contribute to 5km but not at tree level to 9. This can also enhance the predictivity of 
models by reducing the number of parameters, since one can assume that all parameters in 
M°, Mp, M v , M^r, and the right-handed Majorana matrix Mr are real, and that the only 
phase (and only one is needed) comes from A. This is the assumption we shall make in the 
illustrative model we present below. 

Returning to the issue of mass and mixing hierarchies, let us assume that the matrix A 
that characterizes the mixing of 5« with 5^, has the form 

f 1 o o\ 

1 \t 



A 



(7) 



yO \t* A j 



where the real parameter A 1 and the complex parameter t has magnitude of order one. 
As we shall see shortly, it is the smallness of A that gives rise to mb,m T <C m t , while the 
\t\ ~ 1 explains the large atmospheric neutrino mixing. The phase of t, the only phase in 
the model, is what produces the KM phase. We shall see later that the form in Eq. ((7|) is 
easy to obtain. 



To illustrate our basic approach we now present a toy model in which the underlying 
mass matrices have the following simple "textures" : 
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where 5,5' <C t u ,t d <C 1. The similarity of these four matrices is assumed to come from 
5*0(10). In SO{10) one would have the 10; + 5j + lj come from a 16;, whereas the extra 
vector like fermions 5^ + 5- could come from a 10 j. 

The textures in Eq. (JHJ) can be obtained from simple 50 (10) operators. In partic- 
ular, we assume that the 33 elements come from a term of the form ^,33163163(10^). 
Thus, what we have called mu and mo in Eq. (JHJ) are given by mu = h^(H u (10)) , and 
m D — h^(H d (10)) . If the two Higgs doublets of the MSSM came purely from the 10#, i.e. 
if H u = H u (10) and H d = H d (10), then we would have mu/m^ = tan/3. However, one 
expects in a realistic 5O(10) model that H u and H d will come from a mixture of several 
5O(10) Higgs multiplets. Thus, we may write H d = cos7d H d (10) + sm^ d H d (other) and 
H u = cos^u H u (10) + sm.^ u H u (other). Inverting these, we obtain m D = h 33 v cos j d cos (3 
and mu = h^v cos7„ sin/3. Therefore, the usual tan/3 parameter of the MSSM is 
given by tan/? = {mu /m D )(cos^ d / cos7 u ). From Eq. (jSJ) one sees that the top quark 
mass is m t = mu- Therefore, l/cos7 u = h^v/mt) sin/?, and we may write tan/3 = 
(mu /mp) [cos 7 d (v /m t ) sin/?]/i 3 3. The expression in the square brackets is less than or equal 
to 1, and /133, which is a Yukawa coupling in the 5O(10) theory, cannot be much larger than 
1 without destroying the perturbativity of the theory below the Planck scale. Thus, the 
value of the parameter mu/mo, which can be determined by fitting the quark and lepton 
masses, puts an upper bound on tan/?. We shall find that mu/mu, = 2, so with /Z33 = 1.5 
to 2, the value of tan/? is consistent with the experimental lower limit of 3 
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Now, given Eqs. (jUJ), (J7J), and 



one has 
/ 



M. 



D 
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5 5Xt* 
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At A 
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and 
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Mi = 5 3e d At m D . (10) 
5\t* 3e d \t* A j 

Of course, from Eq. (jUJ) one sees that Mjj is already given in Eq. (jSJ) • 

Simply by inspecting these matrices one can observe several significant facts. First, the 
masses of m& and m T are suppressed by the small parameter A, whereas m t is not, so that 
m\>, m T <C m t can be explained without requiring that mu/rriD be extremely large. Second, 
the masses of the first family will be almost unaffected by the parameter A, so that m d 
and m e will not be similarly suppressed compared to m u . Indeed for my/mu of order one, 
m u ~ m d, as observed. Third, there emerges naturally the "lopsided" structure discussed 
in many recent papers |fg. That is, we see that the 23 element of Ml is much larger than 
its 32 element, whereas for Mo the opposite is the case. This comes directly from the 
fact that Md = M D A T whereas Ml = AM". This lopsided structure explains why the 
atmospheric neutrino mixing angle (which gets a contribution from (Ml) 23/ {Ml) 33) is of 
order \t\ ~ 1, whereas the corresponding quark mixing V c t (which gets a contribution from 
(M D ) 23 /(M D )3 3 ) is only of order e d |t| < 1. 

One can read off from the simple forms in Eqs. Q and (jlOj) the following approximate 
relations that hold at the GUT scale: 



m t = mu, 



m u = (S 2 /e u ) mu, 



m b = Aa/1 + \t\ 2 m D , m s = (e d / y/l + \t\ 2 ) m D , m d = (5 2 / e d ) m D , 



m T S X^l+W m D , m M = (3e d /^l + \t\ 2 ) m D , m e = (5 2 /3e d ) m D 
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and 



V cb 



teA) (4) 
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(12) 
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-6' + V us V cb . 



From the form of the Ml (Eq. (|TU)l) and (Eqs. (jSJ) and (JHJl) it can be seen that the 
tangent of the atmospheric neutrino mixing angle is controlled by t, which therefore must 
be of order one. That in turn implies, through the equation for V cb , that V cb ~ m s /m b . 
This succesful qualitative relation between the atmospheric neutrino angle, V cb and m s /m b 
is characteristic of lopsided models. 

In this model there are eight parameters (mu/rriD, e«, 5, 5', A, \t\, and 6 t ) to fit 
twelve quantities (eight mass ratios of charged leptons and quarks, three CKM angles, and 
the KM phase). There are therefore four quantitative predictions, which can be taken to 
be m b = m T , m s = m d = 3m e , and the value of the Cabibbo angle. In addition, 

the atmospheric angle is predicted to be of order one, though it cannot be predicted more 
precisely than that without knowing Mr. 

We can easily determine the approximate values of most of the parameters of the model 
from Eqs. (jllj) and ()12|). We take the values of the quark and lepton masses and CKM 
mixings at the GUT scale to be m t = 112 GeV, m b = 0.96 GeV, m T = 1.16 GeV, m c = 0.27 
GeV, m s = 0.015 GeV, m M = 0.069 GeV, m u = 0.57 MeV, m d = 0.86 MeV, m e = 0.334 
MeV, \V cb \ = 0.0357, and \V US \ = 0.222. These are found by extrapolating experimentally 
determined central values at low scale [9j to the GUT scale using the following procedure. 
First, we propagate the masses of light quarks and leptons from 2 GeV scale to Mz scale 
using the 3-loop QCD and 1-loop QED renormalization group equations (RGEs). Then, we 
perform additional running from Mz to m t scale using the Standard Model RGEs. (The 
relevant renormalization-group (3 functions are summarized in Ref. ^(|.) Finally, assuming 
all SUSY particle masses to be degenerate at m t we run the masses and mixings to the GUT 
scale Mqut ~ 2 x 10 16 GeV using the 2-loop MSSM (3 functions summarized in Ref. jll|. 
In the final running we set tan/? = 3. 

The equation for V cb tells us immediately that \t\ = V cb m b /m s w 2. From Eq. (|TT| one 
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has that (m u m c )/(mdm s ) = [vtiu /mo) 2 a/1 + \A 2 i which implies that mu/mn m 2. The 
equation A = (mb/m t )(mu /mo)/ a/1 + \A 2 then gives A ~ 10~ 2 . 

The value of \e u \ is given approximately by \Jm c /m t = 2.4 x 10~ 3 . The equation = 
|(m At /m T )A(l + |t| 2 ) gives je^l ~ 10 -3 . It is gratifying that e u and come out to be of the 
same order. If we choose the relative sign of e u and to be negative, then we get a good fit 
to the Cabibbo angle: V us ^ y/m d /m 8 {l + |t| 2 )~ 1/4 + ^m u /m c = (0.2) (0.7) + (0.05) = 0.2. 

The value of S is determined from 5 2 = m u m c /m 2 to be 10~ 4 . Finally, the parameter 5' 
and the phase of t can be determined from the real and imaginary parts of V ub . Specifically, 
one has 

V ub /(V us V cb ) = 1 - 5'e^/\V us V cb \. (13) 

One gets a fairly reasonable fit from the following values of the parameters of the model: 
mu /m D = 2.03, A = 1.03 x 10~ 2 , |t| = 1.45, e v = 2.38 x 10~ 3 , e D = -2.14 x 10~ 3 , 
5 = 1.12 x 10~ 4 . The resulting masses and mixings and the experimental values extrapolated 
to the GUT scale are compared in Table |lj It is apparent that the fit is not completely 
satisfactory. In particular, the mass of the r comes out about 15% too small. This is typical 
of grand unified theories. Simple GUTs generally predict m b = m T at the GUT scale, 
whereas the data tend to give m T about 15 to 20% larger than m b . There are a number of 
ways of improving the agreement, including supposing that m b gets corrections from sparticle 
loops. Also off considerably here is m s . The Georgi-Jarlskog relation m s = m T /3 is built 
into a choice of Clebsch in this toy model. But that relation is known to give a value of m s 
that is somewhat large compared to the values favored by recent lattice calculations [l^ . 

While this model does not give a perfect fit, it is simple enough to illustrate the basic idea 
we are proposing in a transparent way. It seems likely that models based on these ideas can 
be found that give better fits. Another possibility that might be realized in this approach 
is a "doubly lopsided" model One could imagine, for example, that the matrix A had 
the form 



I 1 A*i^ 



A 



(14) 



1 Xt 2 
y Xti At<2 A j 

with A < 1 and ~ 1. If the underlying matrix M° has a hierarchical form, with the 
33 element being the largest, then the effective low-energy mass matrix = AM" would 
have the doubly lopsided form, with the 13, 23, and 33 elements all being of the same order. 
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TABLE I: The values of the quark and charged lepton masses and the CKM angles and V us at 
the GUT scale in the model (with parameter values given in text), compared to the experimental 



values extrapolated to the GUT scale. 
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This is known to be able to give in a simple way the correct "bi-large" pattern of neutrino 
mixing angles, with U e % being small 1^ . 

We now turn to the question of whether the form of A given in Eq. can arise naturally. 
Consider the special case where B is diagonal and where the only non-zero elements of A 
are the diagonal elements and the 23 element: = bi5ij Aij = a^ij + 04^2^3- Then it 
is easily found that for b\ < ai, 62 < 02, and 63 > 03, the matrix A has the form given in 
Eq. ((7|) with A = |a 3 | 2 /(|a 3 | 2 + I&3I 2 ), and t = — (b^a^) / (a2\a 3 \ 2 ) . Of course, there are other 
forms of A and B that also give Eq. (|7J). Another simple example is that A is diagonal and 
B has nonzero diagonal elements and 23 element. 

In conclusion, we have found a framework that differs from most "texture" models of 
quark and lepton masses in several respects. First, it can partially explain the fact, usually 
treated as an accident, that m u ~ md,m e , while also giving m t ^> m^,,m M . This it does, not 
by requiring tan f3 to be large, which might be somewhat unnatural, but by mixing the b and 
t strongly with vectorlike fermions at the GUT scale. Second, it combines predictive textures 
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with a structure that realizes a non-axion solution to the strong CP problem proposed many 

n 

years ago |4j. By allowing most of the parameters to be real, even though CP is violated, 
it has the potential of giving very predictive models. And it gives rise naturally to the 
"lopsided" kind of structure that has been proposed to explain the largeness of relative 

to V cb y. 

The toy model we have described illustrates the essential ideas in a transparent way. 
However, it would be good to find a model which is more predictive and which does a better 
job fitting certain quantities, especially m s . It would also be interesting to investigate further 
models of this type that are "doubly lopsided" P- Il^jl. 
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